The code and figures in this folder is for Case 3S (Basin 3 under steady state condition) in Fig. 1
